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A'b8treu:t 


laqpanBions  or   "addition  theorMis"  for  the  spherical  wave  fttnc- 
tlons  J^(ka)  P^»^(co8©)  e^^  .  h^^^(kE)  P^'"^    (coaO)  e^"^  .   and 

h^^^(kR)  p|^''(co80)  •      ^     ,  with  reference  to  the  origin  0,   have  "been  ob- 
tained in  terns  of  spherical  wave  functions  with  reference  to  the  origin 
0*,  iriiere  0*  has  the  coordisates    (r   ,0  , 0     )  with  respect  to  0. 

Using  the  ahOTe- mentioned  addition  theorems,  we  hare  ohtalned 

an  expansion  for     — ,   a  ware  irtilch  is  spherically  symmetric  ahout 

'l2 

the  soaroe  Q,  irtiere  ^  is  referred  to  the  origin  0,  in  terms  of  products 

of  spherical  wares  with  r«ipeet  to  origin  0  and  si^erloal  wares  with 
resp«ct  to  0*. 


-  2  - 


1.  Introduction  ' 

When  a  plane  wave  is  incident  on  a  configuration  of  apheres,  the 
scattered  field  may  "be  ohtained  "by  a  series  of  successive  approximations. 
The  zero-th  order  approximation  i3  just  the  sum  of  the  fields  scattered  hy 
each  individual  sphere  when  the  excitation  of  each  sphere  is  taken  to  he 
the  origineil  plane  wave.  Higher  order  approximations  take  into  account 
contrihutions  to  the  excitation  of  a  particular  sphere  from  the  waves 
which  the  remaining  spheres  scatter.   In  particular,  the  first  order 
approximation  is  the  field  scattered  "by  the  configuration  when  the  ex- 
citation field  is  taJcen  to  he  the  initial  plane  wave  plus  the  zero-th  order 
approximation  to  the  scattered  wave. 

In  order  to  carry  out  this  approximation  it  is  necessary  to  de- 
compose the  scattered  field  of  a  sphere  with  center  at  0  into  incoming 
spherical  waves  for  a  second  sphere  with  center  at  0'.   Such  a  decomposi- 
tion requires  an  "addition  theorem"  which  expresses  a  spherical  wave  with 
center  at  0,  in  terms  of  spherical  waves  with  center  at  0'. 

Such  an  expansion  or  "addition  theorem"  for  cylindrical  waves  is 

1  2 

well  known  .  Using  this  expansion,  Twersky   was  atle  to  calculate  the 

scattered  field  ottained  ty  a  plane  wave  striking  a  configuration  of 

cylinders,  taking  into  account  the  contri"butions  to  the  excitations  of  a 

particular  element  hy  the  radiation  scattered  "by  the  remaining  elements. 

The  analogous  case  of  a  plane  wave  striking  a  configuration  of  spheres  may 

he  treated  with  the  help  of  the  addition  theorems  we  have  ohtained  in  this 

paper. 

The  first  expansion  treated  in  this  paper  is  that  for  the  standing 

spherical  wave  j  (kR)  p'   (cosO)  e  ^      with  reference  to  origin  0.   This 

spherical  wave  is  represented  as  an  integral  of  plane  waves  over  all  possihle 
directions.  These  plane  waves  are  also  referred  to  the  origin  0.  A  trans- 
formation is  then  made  to  ohtain  the  plane  waves  in  terms  of  the  origin  0'. 
In  order  to  evaluate  the  resulting  integral  it  is  necessary  to  have  a  for- 
mula which  expresses  the  product  of  two  associated  Legendre  functions  in 
terms  of  a  sum  of  associated  Legendre  functions.  Once  this  is  ohtained 
hy  using  a  formula  due  to  Infeld  and  Hull,  the  evaluation  of  the  integral 
is  possihle  £ind  the  final  form  of  the  expansion  is  a  sum  of  spherical  waves 
with  reference  to  the  origin  0', 
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war*  _ 


Tha  ■•oond  ezpantlon  treated  la  that  for  the  outgoiog  spherical 
h  '(kH)  p|"^(cotO)  •  ^  vith  reference  to  the  origin  0.  The  same  pro- 


cedure Is  used  as  for  the  expansion  of  j  (kE)  p'°*(cob  O)  e^  .  However 

n      n  * 

in   this  case,  the  integral  is  over  plane  waves  with  real  and  complex  direc- 
tions. This  lead*  to  convergence  difficulties  when  the  transformation  from 
the  origin  0  to  the  origin  0'  is  effected.  These  are  discussed  in  Appendix 
III.   The  final  form  of  the  expansion  as  a  sum  of  spherical  waves  with  refer- 
ence to  the  origin  0'  is  then  ohtalned. 

From  the  first  two  expansions  we  then  ohtaln  an  expansion  for 

h^  (kR)  P^  '(cosO)  e   ^  with  reference  to  origin  0,  in  terms  of  spherical 


wave 


6  referred  to  origin  0'  "by  using  the  relationship  h  HkE)  +  h^^^(kE)=2j  (kfi). 

The  final  prohlem  considered  is  the  expansion  of  the  spherical  wave 

ikr^2 

e 

— ,  where  r,  _  is  the  distance  "between  a  point  of  ohservation  P(r,  ,9,    i,  ) 

•r^2  At  1"    1  "^  1 

and  the  source  point  of  the  spherical  wave  Q,(H,0,  fi  ),   the  coordinates  of  "both 
points  helng  taken  with  reference  to   the  origin  0.      Then, 

— cem  "be  expressed  in  terms  of  a  sum  of  spherical  waves  with  reference 

^2 

to  0.  However,  using  the  results  of  the  precedli^  sections  an  expansion  of 

/^12 

— In  terms  of  products  of  sphericsQ.  wares  with  respeet  td  both  origin  0- 

and  to  origin  0'  is  then  ohtalned.  This  expansion  is  useful  in  the  variational 
treatment  of  diffraction  pro"blems. 

We  wish  to  express  our  thanks  to  Dr.  Victor  Twersky  who  suggested 
this  pro"blem. 
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2.   The  Ejcpansion  for  j    (kR)  p'°'(cob  O)  e^^ 

■^ ~n n 

Consider  a  point  P,  %ftiich  has  the  spherical  coordinates   (R,0,  (5  ) 
vlth  respect  to  the  origin  0.     A  plane  wave  coming  in  along  the  direction 
0  =a  ,0=  3  can  he  expressed  as  follows,   in  terms  of  spherical  waves  with 

center  at  0^^  : 


iia  cosV 


00 


e 


£    i^(2n+l)  jj^(Ml)  P^CcosV),         or  using  the  addition 
n=0 


theorem  for  Legendre  polynomials, 


^IkR  COSV  ^  g      f-         in(2^,^^j^(^3   U-l^  j,|n>|(^^^^)pU  ^^^ 


00        n 

n=0     m=  -n 

Here  V  is  the  angle  hetween  the  directions (  a,  0  )  and  (O,  P  )  so  that 

cosV  =  cosO  cosa  +  sinO  sina  co8()!L^). 

Jrom  the  expansion  (l),   an  integral  representation  of  elementary  spherical  wave 
functions  can  he   ohtained.     Multiplying  hoth  sides  of  the  equation  hy 

p'™'(co«a)e         sina  ,   integrating  over  a,  and  P    ,   and  using  the  orthogonality 

properties  of  the  Legendre  functions,  we  ohtain  the  well-known  formula^ 

4    c,r.\«Ul    f  /^^  ini0-     i~^   /*    /    i^  cosY  |in|/  .    imP     .        .      .- 

J   (kR)P'    '    (cos  0)e     '^=     11^/    /  e  P     '  Uosa)e         sina  da  dp  , 

"'O  *'o     O 

Introdtice  a  new  origin  0',  where  0'  hss  coordinates  (r  ,  6  ,  i) 

with  respect  to  0  (Fig.  1)  and  let  (r,0',j6')  he  the  spherical  coordinates 
of  the  point  P  with  respect  to  0'. 


(2) 


7igur«  1 
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¥e  shall  now  o"btaln  an  aocpansion  for  a  standing  BphericsLL  wave     eo'ctuid.  0 
such  as  Jj^(fcB)  P^"     (cosO)  a     "in  terag  of  standing  spherical  waves  around 

0'     such  as  Jj^(kr)  pj^^'CcosO')  e^^^'. 

Jirst  we  introduce  a  set  of  rectangular  coordinates  around  0  and 
a  parallel  set  around  0'.     Let   (x,y,z)  "be  the  coordinates  of  P  relative  to 
0,  and  let   (x'.y'.z')  he  the  coordinates  of  P  relative  to  0'.     We  have 

«  =  R  cos  0  X  =  R  sin  0  cos   0  y  =  H  sin  0  sin  0 

z»  =  r   cos  0'         x'  =  r   sin  0'   cos   f>  y»   =  r  sin  0'    sin  0 

z  =  z'  +  z^  a  =  I'    +1^  y  =  y,  ^.y^ 

How  since     cosv  =  cos  Q  cos  a  +  sin  0  sina  cos    O-0),   it  follows  that 
R  cosv    =  z  cosa+  x  sinacosP  ■>■  y  sino-sinP 

«   (z»+z^)  cosa^-  (x'-ht^  )  sinacosp  +  {^^  "^  7   )   slnasinP 

=  rQsosO'coaa*  sin  0«   cos^'sinctcos  P+  sin  O'sin^'sin  osin  p^ 
+  r^[co8  O^cosa+  sin  OqCos^^   sinacos  P  +  sin©     sin0     sinOsinP] 


which  can  he  written  as 


R  cosY  =  r  cosY'  +  r     cosY       .  (t) 


Here     v'   is  the  angle  hetween  the  direction  (  a,P  )  and  (O',0'),  while     V^ 
is  the  angle  hetween  the  direction  (a  ,P   )  smd   (0  ,0     ). 

Using  equation  (3)   in  equation   (2)  we  have 

J    (kR)pl'^l(cos   0)e^   =  C  /"/e^l^<^°»^'/^o'°*V|inL  ,     imp  ,      , 

•'n^^-'n     ^°°^      '^  IttTy    J  ®  n   ^cos'^)  ®     ^ainododp.      (h) 

For  the  term  e  we  use  the  expansion  corresponding  to  equation  (1)  : 

IkrcoB'Y'        ^      ^       n, 
e 


=  21      E     i"(2^*l)     fe^4^Jjkr)pU(cosQ.)4^Usa)e-i^^'-P) 
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It  can  be  shown  (see  Appendix  I)  that  the  ahove  expansion  is  a  tmiformly  conver- 
gent series  and  hence  when  it  is  substittited  into  eqtiation  (*»)  the  order  of  suu- 
matlon  and  integration  may  he  interchanged.  Therefore  we  hare 


J^CkH)?!"!  (cos«)e^°^  '^t    t  Ji^a^l){^|4|  J„(lcr)pl^l  (cose-)   (5) 

-iu«Ji  aiT  ,TT  (  ikr.cosY.  alfcl  /    \»U|/    \  i(a«H».)P  .      1  ^  jj 
«  ^^     A  /   J®         ^  '  (cosaJP ' '^ '  (co8a)e    '^  sinal  dadi. 

In  order  to  siniplify  the  above  equation,  it  is  useful  to  hare  a  f omnia 
expressing  the  product  of  two  associated  Legendre  functions  in  terms  of  a  sum  of 
associated  Legendre  functions.  This  foroula  is  obtained  in  the  following  oanner: 
We  first  assume  that 


pl»l(cosa)pi^l(coea)  -  L  ^'T  1°'  ^i""'*'"' 
n  '    '  V     "■         '   e-  P,n,p    p 

where  p  «  i>fn,  v*-n-2,- »'V-n,   the  P'°'(co8a),  p'^'(cosa),  p'**'  '"'are  noroalited 
legendre  functions,  and  eosa  may  take  on  complex  valaes.  Then,  multiplying  the 
above  expansion  by  p'^'  '°' (cosa)Bina  and  integrating  with  respect  to  a,  we  have, 
using  the  orthogonality  properties  of  Legendre  functions 

r   pUKUl  (co8a)pl»l  (coBa)pi^l  (C08a)8inada  -  ^\''\*\^\ . 
i/  0  q.  n        o  q 

However  from  a  paper  bjr  Inf eld  and  Hull'-  ■'  we  have  the  following  formula  ; 


/:pM*1-I  (cos.)f  Ul  (cosa)f  1^1  (cosa)dc.  [tT^]^) Ui^^^l^^^^^ 
r  (g-KluUUI  )y(q-lul-UI  )/(.v|ul)/(n-|m| )/-]  v^ 

ilLtl-UIr  S:^+   I„|+i  (n^|mki).'(.>fq-|m|-i)/ 1 

^       Sr  I  (q-|^|-lm|-i).'tiWn:U|-iS/(.-q*|mki).>j] 


HOW.   pl^l*I»l   .  ll^^l*l»lpl^l*l»l.where  hI^I*!"'-  li^iini  [M±k\\n 


^i/2n._i„_i„iN/-ii/2 


Using  this  relationship  and  similar  ones  for  P'°'  and  P'^'  we  arrive  at  the  formula 
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pl"l  (co8a)pl^^l  (coia)  -  ^  al-"!*'^'  pl»l*l^^l  (cosa)  (6) 

n      '   V  '—    p.n.o    p  ^   ' 

where  p  =  i>*-n,  i>»-n-2,  iH>n-A, ,  u^n,  and  where,  for  any  coefficient  a'°'  l*^'  . 

q,n,p   • 

we  have  the  relationship   i  i . ■  ■ 

n       V 
or 


q.n.P  jjlm|„|pil  q,n,» 


Uklul,    (n^-.>-q-l)//(2q-^l)(o-UI-lul)/  ^  y-      {  (-l)~^ '""'^  V|mkll/(i>fq-|m|-iy 
\^r^,y,  (tt^-q-»)/M«><-q-n)  //(q+iHn+l)//   ^  ^^    |(q-|nil-ln|-i).'(iX(n-|in|-iy(i>-q+|in|  +  iy/ 


(7) 
Here    (8>!J  =  s(8-2)  (s-^).. .   2  or  1.    (o)  ,'/  =  (-!)//=  1,  and  TV-n^  q  <  i>fn. 

(In  Appendix   II  aome  of  the  a'    '    '^'   have  been  calculated  explicitly  for  some  small 

q,n,  17 

values  of  n  and  m. ) 


Using  equation   (6)   in  equation   (5)  we  have 

^.l^^l^Wj  (kr)pl^l(c 

(8) 


j,(«,pi°i (c.,.).*-^=  ^^tt  Er<f:ii[rJ^"'4:i:i'^"-'4'''<-«-)»-'-^' 

y=0  Ll,=— I?  p       >i'-i/.        -fij 


?-i 


where  IT   is  an  integral  of  the  form 

^       Jo      ^  0  ^ 

It  is  possible  to  evaluate  these  integrals  directly'-  J  (compare  equation  ^),  and  we  have 

^•^-  UniPj^(kr^)p],^^l*l«'-(cose^)e^^"^^>^o  .  (9) 

Substituting  equation  (9)  into  (8)  we  have  for  the  expansion  of  a  spherical 
wave  with  respect  to  one  origin,  in  terms  of  the  spherical  waves  with  respect  to  another 
origin,  the  formula 

ML   ~  " 


J  (kE)pi»i (co8«)e^»p=  r.  ml 

wee  \i^-V   p 


^  (HlTTTy^  p,n.p 

Jp  (kr^)j^(kr)p],^l*l"l(cos«^)pl^l(cos6')e^("^^^^^oe-^^^'J  •  (10) 
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3.     The  ExpaMlon  of  h^^^lcH)pl'°lco8«)e^'"^  and  h^^^kE)?'"' (0086)6^°^ 
an  n  n 

In  order  to  o^btain  an  addition  theorem  for  the  outgoing-  spherical  wave 

h'  ^(kS)P^"' (co8e)e  "^  we  must  first  get  the  Inte  ral  representation  for  this  fTinctlon. 

We  have  the  formula'-  J 

"^  "^    Jo    K         e^^l'^^^^l  slnadad^  (11) 

where  r. ,  «,,  fi.   are  the  coordinates  of  a  point  P.  with  respect  to  the  origin  0,.  Now, 
let  the  point  P^  have  spherical  coordinates  (B,6,^)  with  respect  to  the  origin  0,  and 
let  0.  have  coordinates  r  ,  ^  ,  ©  with  respect  to  0«  From  equation  (3)  we  see  that 

BcosV  =  r°cosY°  +  r.cosVj^ 

or  ikr^l^cosY^  _  -lkr°cosV°  IkEcosV  (12) 

where  as  "before  cosV  =  cosacos6  +  slnctBln€co8(3-0). 
Now  from  equation  (l)  we  have 


-ikr°cosV° 


n=o  m=-n 

Therefore,  using  the  formula^  -^ 
ikri   00  n 


=  f  i:  {r^2n.l)  4^W4  J  (kr°)pl»l(cos«°)pl»l(cosa)e^"'(M']: 
n=o  m=-n  (I3) 


1         n=l  iii=-n 
and  equation  (12)  and  (13)  in  eq-uation  (11)  and  comparing  the  coefficients  of  «acu 
j  (kr  )p''°'(cose  )  we  find  the  Integral  representation  to  be   '-  ■' 

Tt  ,  ^^ 

v(l)/,,T5^•DUI  c  a\   M  i"'^  r        f^       IkEcosV-UL    a  im^  ,   ^  .«        f^c\ 
h   '(kR)P'  '  (co8e)e  ^=   — r-  J  e  P'  '  (co8a)e   '^ainadodP.        (15) 

n       n  iiTT  "  0   >'  o  n 

Now  let  the  point  P  have  spherical  coordinates  (r',9',0')  with  reference  to 
another  origin  0',  where  0*  has  coordinates  (r  ,6  ,0  )  with  respect  to  0.   (See  Pig.  1), 
Prom  equations  (l)  and  (3)  we  see  that 

Rcos"V  =  rcosV '  +  r  cosV 

o    o 


and  thft 

IkrcosY' 


e 

0=0  \i,=-0 


r  E  i^'Ca^DJ^Wjf  J^(kr)pN(co8e.)pW(cosa)e-^^(^'-P^ 


Whan  these  two  relations  ^r®  substituted  Into  equation  (15) »   it  can  he  shown  (see 
Appendix  III)  that  It  is  possible  to  exchange  the  order  of  summation  and  Integration 


provided  that  r  <t   .     Therefore  we  have 

0 


a 


(kE)pl°l(co8e)e^"'^-^E  t    ll''(2.>fl)fe|4^J^(kr)P^^I(cos6')e-^^^' 

0=0  (j,=-»[  ' 

J^jl  -'>r,ce8e,p|^|(^^3^jp|ml  ^^^^j^K^^m)?^^      ]  ^   ^6) 


for  r  <  r. 
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Fov  we  have  the  fonrala  (see  equation  (6)) 

where  p  «  »»-n,  W-n-2,  •••  v-n  euad  the  a'^'  '  *  are  glTen  hy  equation  (7).   Subetituting 

P»D.|  P 

the  a¥oTe  esiMiattlon  Into  equation  (1?)  we  hare 


h^^^(kR)p''"'(co8«) 
n      n 


ia.n/x'     vcuB'v/o    — — >!—  /    /      /J*  K^-""^-*-  J~f — n — I  \  J  ■»     _   J 

4 . .  /  I  , .   _ 

r<r^,     (18) 


where  the  xj^'"  are  integrals  of  the  form 

^•°  -  /f/o  "  "^  e^'^o°o««op|;^I^M  (cosa)e^^^*-^^inadadP. 
These  integrals  can  be  evaluated  directly.  Using  equation  (15)  we  find  that 

kJ'"=  2niPh^^^(kr^)p],^*l*l'"le^^"^^^o- 


Therefore,  equation  (l8)  becomes 

00   V 


pl^^l-'l»l(cose  )pl^l(co8e')e^^"^^^o"^^i 


,_,._..  ,         r  <  r 

p      '    0'  »  J  0 

When  r  <  r.the  resulting  expansion  would  he  the  above  equation  with  r  amd  r  interchanged, 
Eiert 

(1), 


0 

Therefore  we  have 


h 


u=o  M^-o  pi 

pl^*l*f»l(cose  )pi^^l(cose')e^^"^^^o  e-^^<*'  .  (19) 

p  0   »  J 

where  r^  is  the  saaller  of  r  and  r  and  r^  is  the  larger  of  r  and  r^, 

fhls  is  the  expansion  for  a  spherical  wave  h^  '  (kE)?'"*' (coBe)e  ^  about  one  origin, 

in  terms  of  spherical  waves  around  another  origin. 

To  find  the  expansion  for  h^  '(kE)p''°' (co8«)e  "^  we  use  the  relationship 

n       n 

h^^^(kH)  +  h^^^kH)  =  2J^(kl) 


and  the  expansions  for  J  (kR)p''°' (co80)e^°'^  and  hy^lcB)p'"' (0089)6^""^.  This  gives  us 
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h<^)to)pi»i(=o..).'<E  t  e{^(2-i)>i:!,"|j,(^<)»^^>(^>) 


p  0       1» 


(20) 


ikr 


12 


^.     The  ExpanBion  of 


12 


Let  P(r.  ,6,  ,0, )  be  a  point  of  ©"beervatlon,  and  Q(E,e,^)  "be  the  source  point 
of  a  Bpherical  wave  where  the  coordinates  are  referred  to  a  fixed  coordinate  system 
whose  origin  is  at  0.  [jig.  2}  ,     The  distance  from  Q,  to  P  is  therefore 


^^12  ~  *  ^^'*"'^1  -  2Er^co8V  ^ 
where  cos'v  =  cose  cosO,  +  sinesinQ,cos()d-^. ).  Hence  we  have  the  expansion 

lkr^2    00 

2.^^ =1:  (2iH-l)P^(cosY)j^(kr Jh^^^kr  ), 

12     n=o 

where  r^  is  the  lesser  of  E  and  r, ,  and  r^  is  the  greater  of  E  and  r, , 


Q(R,^,«^) 


l9] 


(21) 


Tig,   2. 
Using  the  following  addition  fornnila  : 

P(cosV)=L      l^rHlf    P'"' (cose )pl»l  (cose, )e-^»(M> 
n  * —       (n+  m   )J       n  n  1 


mF— n 


in  eqtietion   (21)  we  have 
12        00     n 


i-^.t    t     (2n*l)4M44j „(krJh(^'te.)pl»l(co.e)pl"l(co.e,).-^"(»'l-'''- 
r,  „    * —  * —        (n+  m  ;Jn<n     >n        n       x 
12    n=o  m=-n         ' 

Consider  now  the  spherical  coordinates  {1^,9^,^^)   of  the  point  q  with  respect 
to  the  origin  0',  where  0*   has  coordinates  (r  i^q.^q)  wi*^  respect  to  0.  [See  Fig.  2] 
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0'(roA.<f>o) 


yig.  3. 

Then  we  can  apply  to  eqtiation  (20)  the  following  erpanslona  1 

J  (idDpl-l  (coBe)e^<  E  t    n 


v=o  ^i=-»  p 


(wU|)l     -^p^^o^^u^^^Z^^p.n.u  ^p     ^ 
pW(cos92)e^^"^^^Oe-i^^2  1  ,  (23) 


COS 


^CPJ^^I-^I"!  (co9«Jp[^^l(coee2)e^^'^^^o«-i^^^2J;       (2l*) 

where  r'  is  the  lesser  of  r„  and  r  and  rJ  is  the  greater  of  r^  and  r  .   The  final  re- 
<  2      o)     >  2      0 

Bult   ,  after  substituting  either  equation  (23)  or  (2^)  into  (22)^  is 

ilcr 

--r     erf:   i:B'^''l'"l(r,.r  ;r9  ,0  )pl"'l(coe9,)pl^^l(cos0je-^("°^l**^^2^ 

r,     ^—  ' —  * —   * —   ^—  n.tf.p   ^  1  2  0  0  '^o  n       1   u       2 
12    n=o  1^0  in=-n  ^=-»  p     "^  /„<) 

where  p=i>t-n,i>fn-2,  •  •  .y-n,  and  where  for  r,  >  H 

Bl^l'l^^r^.r.jre   .0   ).,'>^P-n(2n^lj(2;>^l}(n-|n>l  )l(->-|^^l  )i  J^KN j    (^   )h(^^]cr,  ) 
1*   2*   o  o  '^o  (n+|m|  JKiX-Im-I  )i  P.n,»     "p       on  1 


n,»,-p 


and  for  r.  <  H 


Xj^(kr2)pl^l*l»l(cose^)e^^»*^^^o 


Bl^l'l°'l(r     r-r©     6   x   .»H-p-n(2n4-l)(2»>flj(n-Uh.' (l>-|^.|  ).'  al^l-^l^lh^^^kr' )J    (kr   ) 
n.p.p      ^''l''^2'Vo*Po^^  (n+|in|)l(i>f|^l)J  ^    P.n,.;       p      ^'"yJn^^^l^ 

O^(kr4.)p'^^l*l"l(cose^)e'^"^^^^o- 


If  we  let 


Al^l'l"l(r,iO   .^   )  ,     ,>.p-n(2n^l)(n-M)l    (2i>fl)  (i>-|^x|  )i  J^xK|m|p|,.K|m|  ^^i(m^^)P, 

n.u.p     ^   Z*   o'^o'  (n+|m|)l  (i>f|n])i  P.n.u       p  ^  o' 


then  for  H  >  r. 
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l^h^^^rirr'  )^    (Vr    u' 


=i:i:r'''r'2"o-«o-'o''j.(^iV  <'->'^n<^i>^i:i:; '<'a^«o-'o'. 


and  for  H  <  r. 


(1), 


»;:;:i '  ^^v^2'^o-\-(o>h^^z'>i/^.K'^^iKyr  ^^r\-v  • 
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Appendix  I, 

To  show  that   the  expansion 


^IkrcosV.,^    ^     li»(2^l)j^|^j^(kr)p],^^l(cos«.)pl^l(cosa)e-^<'*'-P^ 


is  a  uniformly  convergent  serleB^we  first  write  it  in  the  form 
^Ikrcosv.  ^^  i"(2i>fl)j^(kr)P^(cosV«). 

How  we  know  that  f  or  0  <  Y '  <  2tt,  |p„(co8V«)|<  1,  and  also  |e"^^^*'~^^|  =  1, 
Tor  large  values  of  o,  !j  (ter)|  "behaves  like  (-j-)  .  Hence  we  have 


H   i''(2i>fl)j Jkr)P^(co8V)|  ^  ^  (2^l)(^)^ 


1>=0 


»=0 


W 


00 

0ESO 


However,  ^  (2j>*-1)  (-r-)   is  a  convergent  series,  and  therefore  it  follows  that 


00 

JZ  l''(2i>fl)j^(kr)Py(cos'Y«) 

1^0 


is  a  -uniformly  convergent  series. 


Appendix  II. 

For  general  values  of  o,  n.  ^  and  |ml,a  formula  has  heen  obtained  (see 

equation  7)  for  the  coefficients  a'^'   °  .  Using  this  formula,  we  have  calculated 

p,n,P 

the  coefficients  for  some  special  values  of  n  and  m.   It  is  sufficient  to  evaluate 

the  a'^'  '™'  for  positive  values  of  m  only,  since  from  these  the  coefficients  for 
p,n,  1/ 

negative  values  of  m  are  readily  ohtainahle  and  depend  upon  the  definition  of  P~°(i) 
in  terms  of  P^(x). 

We  see  that  for  n=l  the  a'^'  '  '  =  0  unless  p  =  iH-1  or  i>-l  ; 

Pi  If  " 


n=»l 


m=0 


m=l 


jM.l+|mI 
Vl.l,i» 


(i>-Im.I-»-1) 
(2i>fl) 


2i>fl 


i>-i,i,p 


2i>fl 


-1 
2»fl 


Tot  n=2,   the  a'^j,'*''"'  =  0  unless  p=i>f2,w,  v-2 


n=2 


in=0 


m=l 


m=2 


\>*-2,2,v 


2  lH^^liil+2}lH=Mt2J 


2      (2i>f3)(2i>*-l) 


■^   (2i>f3)(2»><-l) 


(2i>f3)(2i>M) 


lul+lml 
tf.2,u 


(21H- 


H-3)(2u-l) 


1  i^M) 

2    (2i>f3)(2i>-l) 


-6 


(2i>f3)(2i>-l) 


''-2,2,u 


When  n=3,   a'^'"*"'°'=  0  tuilesB  p  =   i>f3,    i>fl,    i>-l,    i>-3. 
P » -'' » ^ 


3  (^lul)(i>fluUl) 
2  (2i>fl)(2i>-l) 


-^  (2i>*-l)(2i>-l) 


(2i>*-l)(2i?-l) 


Tor  idfO 


M-^W\     _  i  (i>-|uk3)(>N.|uk2)(.^|^t|-H) 
»>^3,3,w   "  2   (2i>f5)(2i>f3)(2i>4-l) 


.|^xt+|ml   _  2  (»^!ukl)(u^^-2i>^5LI^) 
i>4-l,3,y  ~  2    (2i>4-5)(2i>*-l)(2i>-l) 


i>-l,3,»     2   (2i>f3)(2iM-l)(2tf-3) 


M-^\m\     ^     1  (i>f|uh(i>f|ul-l)(i>*-lul-2) 
"  -       2   (2i>fl)(2i>-l)(2iv3) 


>'-3t3.i' 


Appendix  III . 

In  order  to  prove  the  validity  of  exchanging  the  order  of  summation  and 
integration  in  equation  (l6),  we  first  write  it  in  the  form 


h^^^lcR)p'"'(c08e)e^<  ^   lim  /^  f  ^        { 
n^'n  2tt   ^^ovq 


00 

51  i''(2i>fi)j^^(icH)P^(co8V« 
u=o 


3 


^ikr^cosV^  P^'°l(cosa)e^"^inadadP  . 


(26) 
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iH)'' 


HoweTer,    Ij^Ckr)!   <    -^      and   IPy(c08Y0|  <     |    cosV  + /cos^'-l   |"     ; 


hence 


00 


^ 


kr 


(2i>fl)|    2     (cobV«Vcob^'-1    I*^ 


where  the  right  hand  aide  is  a  convergent  series,  and  therefore 

CD 

H  i''(2iHl)j  (kr)P^(co8V') 

is  a  •uniformly  convergent  series  for  any  finite  values  of  r  and  V. 

Now,  in  the  finite  region  containing  the  segment  of  the  contour  of  integra- 
tion included  "between  and  0  and  ■x  -  iC,   e   o  °  op'"' (co8a)8inae    is  bounded  and 

<c  n  ' 

hence  the  series 

^ikroCosY^  p['"l(c08a)sinae*"P  2Ii^(2.^1)jy(kr)P^(co8V« ) 
is  uniformly  convergent  in  this  region.  Therefore  we  can  write  equation  (26)  as 

h^^\kR)p["'l(co80)e^<i^liB  )  ^l''(2»fl)j^(la-)/^^2"  P^(co8V  )P^°'' (cosa) 

0— >oo  1  1^0 


ikr.coeV^  imP^ .   ,  . a  ( 
e   0    0  e   Blnctdadp  >  * 


Equation  (2'')  can  then  "be  put  in  the  form 

h^^\kR)p]^"'  (co89)e^"^=^  r  i''(2i>fl)j^(kr)^^^'p"^(cosV<)p['°'  (cosa) 


2tt  -T-ioo 


1>=0 

ikToCOsVQ  jjjg 


e       sinododP 


(27) 


.-n  \    00  2Tt  ^r  -i  00  l     I  llrr    COIV       *     a 

i—  lim     iYl^''i2^l)i    (^)f   L  P,(co3V)pl"''(co8a)e       °         %^"'^inadadP 


J-iC 


If  we  can  now  show   that 
.-n  I    00 


(28) 


2Tt   rr-i  CO  II  ikr  cosv 

-2        _    ,       .....Iml,  s  0  Oe^'^^inododP 


K  =  ir^  lim     \Y1  i''(2i>4.l)j^(kr)/  /^        Py(co8v  -  )P^°I  (cosa)e 
0->oo  1  v=o  ^     " 


=0, 
(29) 


-  16  - 

Ihen  the  proof  will  "be  completed. 

To  prove  that  K  — >  0  as   C  — >oo,  we  first  su'bstitute  Into  epilation   (29)   the 
eTpanslon 

P„(cosv.)=i:   j^Hji  .  pI^^I(cose.)pW(co8a)e-^^^^'-^^ 

(Since  this  is  a  finite  s-am,  we  may  interchange  the  order  of  summation  and  integration.) 
This  gives  ufl 


iter  cosY^  . /_^  v- 
e   0    Oe^^"^^^^^inadadP 


1- 


(30) 


Uow,  in  the  (a, 3)  plane  we  perform  a  rotation  snch  that 

cosa'  =  cosY  =co8cLCOBfi  +Binasin6_co8  (P-ji^) 
o  0  0        0 

e-  =  Mo  . 

Then  cosa  =  cosa'cos©  +8ina'8inG  co8(P'+0  )  and  the  limits  of  integration  remain  the 

0  0  0 

same.     Hence  K  'becomes 

C->oo   ^i>=o  p.=-» 


2tt     ■T'-i  oof  ikr  cosa' 


f    f^         e   °    p'^^'(co8a'co»e  +8ina'8ine  cos(P«+^^)) 

J       J  O  O  0  0 


(31) 


"   i-ic 


p''"'(co8a'co8e  +sina'8ine  co8(p«+0  ))e^^"*'^^  ^^'"^o^sina'da'd^ 
n  0        o       "^o  JJ 


Now  for  Isrge  values  of  \%\,  \ie   have  the  following  asymptotic  form  for  P'  '(z)^   "^ 

,  ,     2Pr(i4)u)^ 

^  (p-UDiRl) 

If  we  let  a'  =  ?  -  iy,  where  C  <  vp  <  oo,  then  for  very  large  vj/,  cosa'  hehaves  like 
ie^  md  sina'  "behaves  like  e]^  .  Hence,  using  the  above  asymptotic  form  we  find  that 


2  2 

P'^l (cosa'cose  +8ina'sine  co8(P'+0^)) 

MX  O  O  O  ' 


p'"' ('coS'i'cosQ  +8ina>8ine  C08(P'+©^)') 
n   \        0         0       0  / 


2T(4)(f)' 


(«^ki)jr(|) 

2^r(n4)(^)'' 


(n-!m| )•[(!) 


-  17  - 


Therefore,  on  substituting  the  ahove  expressions  into  equetion  (29),  and  then  taking 
the  absolute  value  of  K  and  integrating  over  p',  we  have 


K  < 


1-^(4) 


2nr(|)r(|)(n-lm|). 


00   17 

1^  {Yin 

C— >oo   u=o  \x=-v 


(^1.1).     J,(kr)pl,^'(cos6')e-^^^ 


"/  (e^) 


(32) 


The  integral  in  equation   (30)   is  now  a  portion  of   the  real  axis.      Letting 


-V 


t»«  the  integral  hecomes 

00 


"e 


-00 


-kr 


iM-n   2 


0  t 


/  (6^)*^%  "^         e^dn,  =  /  t-"e  '   dt 


(33) 


where  D  =  e  and  therefore  as  C— >oo,  IH>oo.   Integrating  equation  (33)  by  parts, 
it  becomes 


-00 


kT 


S  ' 


i>*-n 


£  t 


dt 


2t   e 


kr 
_  o 


kr 


,  2lgnly.-^>.n-V  2  \, 


(3^) 


Therefore 

.00 


y^  (e'^)   e       e'^dv<  ^r ; 


(95) 


and  using  this  inequality  with  equation  (32),  the  result  la 


K  < 


r(4) 


rlNn/l« 


r(f)r(t)(n-!ml)Jkr^  D->oo 


■^  00   00  2(i>fi)r(^^)         I  I 
lim  D^  2   X:  r     ('^|^|)f  |j,(kr)||pl^'(cosQ.)|D^ 


1^0  M^O 


However,  |p'*^' (cose' )  |  <  ^ 


(kl)" 
and   Ij^(kr)  |  ^^ 


for  large  values  of  u. 


Putting  these  values  for  |p'^' (cosB' ) |  and  |j^(kr)|  into  equation  (3^)  and  writing 


('^)r(»^2)  as  n^).  we  i»"ve 


K  < 


r3/2^j.  y2sine'(n-|m|  )i  D->oo 


-krn  _,      ^krPxP 
lim   D%  2    i:   -2 


1^0 


w; 


I   « 


(37) 
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Eowerer,  2_    r~    "  •         ♦  '^^  therefore  eqijatlon   (37)  'becoifteB 

IXKAU.      ."J''-'°*  (38) 

D->flD 

where  A  Is  the  conetant 


If  we  now  let  D  -^cd  we  see  that  \t\-^  0  prorided  that  r  <  r  ,  Therefore  equation  (28) 
hecomeB 


2"  -r^loo 


h^^^kE)pJ^"l(c08«)e^<  ^  r  i^2tH.l)j^(kr)yr  ^2*  "p(co8V  .)p^»l  (cosa) 

l>=0 

Ifcr  coB^  .  fl 
•   **    S^^'slnodadP      when  r  <  r^        (39) 

and  when  P  (coaV)   1b  excanded  in  terse  of  p'^' (coBa)p'^' (cob«')  we  have  eq-'oatlon   (16). 
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